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ABSTRACT 

By combining the newly infrared photometric data from the All-Sky Data Release of the 
Wide-field Infrared Survey Explorer with the spectroscopic data from the Seventh Data Re- 
lease of the Sloan Digital Sky Survey, we study the covering factor of warm dust (CFwd) for a 
large quasar sample, as well as the relations between CFwd and other physical parameters of 
quasars. We find a strong correlation between the flux ratio of mid-infrared to near-ultraviolet 
and the slope of near-ultraviolet spectra, which is interpreted as the dust extinction effect. 
After correcting for the dust extinction utilizing the above correlation, we examine the re- 
lations between CFwd and AGN properties: bolometric luminosity (iboi), black hole mass 
(i\/BH) and Eddington ratio (L/L^dd)- We confirm the anti-correlation between CFwd and 
Lboi- Further we find that CFwd is anti-correlated with but is independent of L/L-^dd- 

Radio-loud quasars are found to follow the same correlations as for radio-quiet quasars. Monte 
Carlo simulations show that the anisotropy of UV-optical continuum of the accretion disc can 
significantly affect, but is less likely to dominate the CFwD-^boi correlation. 
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1 INTRODUCTION 

The terminology 'active galactic nucleus' (AGN) generally refers 
to energetic phenomena in the central regions of galaxies that can 
be attributed to the accretion of gas onto supermassive black holes. 
Due to its angular momentum, the in-falling gas forms a disc 
around the black hole and slowly spirals inwards as friction torque 
transports the angular momentum out. During this process, gas in 
the disc is heated to a few 10^ K and produces strong UV to op- 
tical continuum. Strong UV radiation ionizes gas surrounding the 
black hole, giving rise to prominent broad emission lines within one 
parsec and narrow emission lines further out. These emission com- 
ponents are further modified by dust in a torus outside the broad 
line region (BLR) and in the interstellar medium of galaxies. Dust 
absorbs and scatters UV and optical light and re-emits at infrared 
(IR) bands, altering spectral energy distribution (SED) as well as 
causing anisotropic obscuration. In the past thirty years, it has been 
established that a large portion of observed diversities, such as pres- 
ence or absence of broad lines and prominent non-stellar contin- 
uum, can be attributed to t he anisotropic obscuration of the dusty 
torus (see I Antonucci|[l993l for a review). 

The study of dusty torus is of great significance in our un- 
derstanding of physical processes of AGNs. The torus is a natu- 
ral reservoir of gas supply to the accretion disc, and thus provides 
a possible connection between the accretion disc at sub-parsec 
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scale and the more extended nuclear disc (e .g. Krolik & Begelmanl 
1 19861 : iHopkins et aLll2005l : lOi Matteo et alj|2005l) . It has been ar- 
gued that the torus is a smooth continuation of the BLR and the 
bou ndary between them is determined only by dust sublimation 
(e.g. Elitzui 20081). S o the torus may define the extension of BLR 
jNetzer&Laodl 19931) . Also, dusty torus may be the source of out- 
flows on the scales of parsecs. The anisotropic obscuration of ion- 
izing continuum causes an ionization cone structure in the nar- 
row line region, which has been detected in many obscured nearby 
Seyfert galaxies. 

Due to its great importance, the geometry, column density 
distribution and the composition of dusty torus have been sub- 
jects of intensive study in the past thirty years. The average cov- 
ering fraction of torus could be constrained from the number ra- 
tios of obscured to unobscured AGNs in the local un i verse. For 
nearby Seyfert galaxies and radio galaxies, ILu et al.l (1201 Oh re- 
ported a Type 2 to Type 1 ratio around 2:1 to 3:1. This suggests 
that an average opening angle of the dust torus is about 45 degree, 
which is basi cally consisten t with that derived from quasi- steller 
radio sources (lB arthelll989h and d irect imaging o f ionizing cone in 
many nearby Seyfert ga laxies (e.g. lWilson & Ulv e stad 1983; Pogge 
' l98^:ISako etal.ll2000l) . While in the literature. [Lawrence & Elvis 

1:1, and so did iReves et al 



20101) suggested this ratio to be 

20081) by reporting approximately equal space densities of ob 



scured and unobscured quasars. Yet a consensus has still to be 
reached on whether and how the open ing angle depends on other 
properties of AGNs. lLawrenc3 (Il99lh suggested that the fraction 



© 2013 RAS 



2 Ma& Wang 



of type I AGNs increases with the increa se of the nuclea r lu- 
minosity. Simi l ar res ults were obtained by ISimpsonl dlOOSl) and 
iHao & Strausl ( |2004 for a large sample of AGNs d erived from 
Sloan Digital Sky Survey (S PSS), and bv iHasingeJ feOOSl) from 
deep X-ray surveys. However, ILu et al.l (1201 ol) did not find such a 
congelation for radio loud AG Ns by taking into account of various 
selection effects, and nor did [Lawrence & ElvisI ( |2010|) for radio 
quiet AGNs. This inconsistency is largely due to the correction of 
complicated selection effects in dealing with the AGNs in each lu- 
minosity bin. 

The infrared emission is a powerful probe of the dusty torus. 
With a maximum temperature around ~ 2000 K, the emission of 
dusty torus is mainly at IR bands. The total IR luminosity from 
the torus is a measurement of the amount of optical to UV light 
absorbed by the torus, while individual features are used to probe 
the composition of grains and their distribution. The Spitzer Space 
Telescope and more recent Herschel provide very rich data-sets, 
which allow us to study different features and parameters of the 
dusty torus. Meanwhile, large unbiased mid-infrared surveys allow 
to study the statistical properties of the dust torus. One of the prob- 
lems that we are primarily concerned is the dust covering factor 
(CF) and its correlations with AGN properties, such as bolometric 
luminosity, black hole mass and accretion rates. In the literature, 
several authors have suggested that the CF d ecreases with increas- 
ing bolometric lum inosity (Lboi) (see, e.g. iMaiolino et al]|2007l ; 
lTreisteretaLll2008l) . However, these studies used either a relative 
small sample, or a combination of a low-redshift, low-luminosity 
sample with a high-redshift, high-luminosity sample. In the latter 
case evolution effect cannot be excluded. 



Recently, the Wide-field Infrared Survey Explorer (WISE; 
IWright et al.ll201(3l) provided a large amount of photometric data 
at near- and mid-infrared bands that can be used to study the 
IR emission as well as dust CF of AGNs for a quite large sam- 
ple. For example, recen t works of iMor & TrakhtenbrotI (120111) and 
ICalderone et al.l ( 1201 2h reported that CF anti-correlates with Lboi 
by studying large magnitude-limited samples. In this paper, we will 
utilize the most recent WISE All-Sky Data Release, combined with 
the spectroscopic data from the S eventh Data Release of the Sloan 
Digital Sky Survey (SDSS/DR7; lAbazaiian etal]|2009l) , to study 
the warm dust emission of 16275 quasars in a redshift range of 
0.76 < z < 1.17 (including ~ 10% radio loud quasars). The sam- 
ple is described in §|2T| We fit the SDSS spectra with power-law 
function to derive the flux in the near-ultraviolet (NUV) range of 
rest frame wavelength from 2000 A to 4000 A (Fnuv ). This NUV 
band is close to the pea k SED of the big blue bump jZheng et al.l 
ll997l : IShang etal.l201ll) . We also calculate the flux in near- to mid- 
infrared (MIR) range of rest frame wavelength from 3 /im to 10 /im 
(^MiR, ^ |2.2| l. which is dominated by warm, and a small portion of 
hot, dust. We find a strong correlation between -Fmir/Jnuv and 
NUV spectral slope a, which can be best explained as dust ex- 
tinction and reddening in NUV (§[3](. After correcting for the dust 
extinction, we estimate the covering factor of warm dust CFwd, 
and examine its relations with other AGN properties (§|4l(. In the 
end, we compare our results with previous works and investigate 
how anisotropic continuum can affect the observed CFwD-^boi 
correlation (§[5](. 



In this paper, we use a flat cosmology with Ho = 
70kms"^ Mpc"\ = 0.3 and =0.7. 



2 THE SAMPLE AND MEASUREMENT OF FLUX 
2.1 Sample Description 

We construct our q uasar sample by cros s-conelating the SDSS 
DR7 quasar catalog dSchneider et al .120101) with the WISE All-Sky 
Data Release catalog. This leads to a sample of 101639 quasars 
with reliable WISE detections in at least one band (96% of whole 
SDSS DR7 quasar catalog). To get a more reliable estimate of the 
bolometric luminosity and mid-infrared luminosity, we require that 
the SDSS spectra cover the rest-frame wavelengths from 2000- 
4000 AA and that the WISE bands access to at least lO/xm in the 
quasar rest frame. This leads to a redshift cut 0.76 ^ z ^ 1.17. 
17639 quasars fall in this redshift range. 

We fit the rest-frame SDSS spect ra, c orrected for the G alac- 
tic extinction using th e dust map of ISchl egel et al., i ll998l) and 
the extinction curve of iFitzpatrickl ( Il999l) , in three relative line- 
free windows 2150-220 0, 3030-3100 & 4150-4250 AA (see 
IVanden Berk et alj200lh with a power-law function 



Fx = c (A/3000)", 



(1) 



where normalization factor c and spectral slope a are two free pa- 
rameters to be determined by standard minimization approaclQ. 
Bad pixels are removed using the mask bits in the SDSS spec- 
trum. To achieve reliable fitting, we further require that the num- 
ber of good pixels must be greater than 50, and that the median 
signal-to-noise ratio (S/N) should be greater than 3 in each win- 
dow. 1220 sources (7%) fail to meet these criteria and are thus 
dropp ed. Another 71 so urces with broad absorption lines (accord- 
ing to lShen et al.ll20I ll) were also excluded. Furthermore, we dis- 
count another 69 sources with fitted spectral slope a greater than 
because their continua are severely reddened and are likely con- 
taminated by starli ght. By this stage, 16275 sources (92%) are left. 

According to lShen et aO ( 1201 ll) , 15279 of these sources lie 
in the footprint of th e Faint Images of the Radio Sky at Twenty- 
Centimeters (FIRST; Iwiiite et al.ll 1 99?!) survey. And 1568 of them 
are classified as 'radio-loud quasars'. In the following analysis, 
we will examine the potential differences between radio-loud and 
radio-quiet quasars on the properties we are discussing. Further- 
more, since our purpose is to study the mid-infrared emission of 
warm dust, we will first restrict our analysis to the primary sample 
of 12483 out of 16275 quasars that have reliable detections in ALL 
four WISE bands, and then examine the properties of the rest 3792 
quasars which only have upper limits of W3 or/and W4 magni- 
tudes. All information for the entire sample relevant to this paper is 
available online as the supplementary material (see Appendix and 
Table[T]for a description). 

2.2 Near-UV and Mid-Infrared flux 

We calculate the integrated NUV continuum flux, -Fnuv, in the 
rest-frame wavelength range of 2000-4000 AA by simply inte- 
grating the best fitted power-law model. It excludes the reprocessed 
components such as B aimer and pseudo Fe II continua, and broad 
emission lines. The enors of -Fnuv come from two sources: the 
uncertainties of fitting parameters c and a, and the uncertainty of 
SDSS spectrophotometric calibration. Using enor propagation, we 
estimate that the former introduces only a small enor, typically 
0.5%. We expect the latter to be the dominant source. To quantify 
this, we notice that the difference between spectrophotometric and 



We use the MPFIT package for all nonlinear fitting in this work. 
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photometric magnitudes has a la scatter of 0.05 mag at r bancQ, 
which is mostly attributed to the uncertainty in spectroscopic cali- 
bration. This indicates the spectral flux at r band has a relative un- 
certainty of ~ 0.047. We will quote 5% in -Fnuv as the calibration 
uncertainty hereafter. 

For mid-infrared (MIR) emission, we will focus on the in- 
tegrated flux in the quasar's rest-frame wavelength range of 3- 
10 jitm. This regime is covered by the WISE W2, W3 & W4 bands 
for our chosen redshift range. Furthermore, the infrared emission 
in this regime is primarily dominated by the warm and hot dust 
heated by quasars. Finally, the mean qua sar SEP in this wavele ngth 
range is relatively smooth according to [Richards et al. liooi) and 
IShangetaLlJioT ij). Therefore, it is relatively straightforward to ob- 
tain the integrated flux from the WISE photometric data without 
invoking detailed dust modelling. 

In practice, we approximate the MIR SED as a broken power- 
law with broken wavelength at the effective wavelength of W3 band 
at the rest frame of the quasar. Then we integrate the model from 
rest-frame 3 /im to 10 /im, to obtain the integrated MIR flux -Fmir. 
To estimate the uncertainty of this measured MIR fl ux, we apply 
the sa me approach to every object in the sample of IShang et al.l 
( l201lh . which have the Spitzer MIR spectra. We find that our ap- 
proach gives -Fmir 3% higher on average than the direct integration 
of IR spectra over the coiTespondent wavelength range, with a scat- 
ter of 8%. Hereafter, we adopt 3% as the global offset and 8% as 
the uncertainty introduced by our method. Finally, we combine this 
uncertainty with statistical uncertainty comes from the uncertainty 
of WISE magnitudes as the final uncertainty of Fmir. The spectra 
slope Q, NUV flux -Fnuv and MIR flux -Fmir for the entire sample 
can be found in the online table (see Table[T]for a reference). 



3 EXTINCTION CORRECTION 

3.1 Correlation between Fmir/Fnuv and Spectral Slope a 

The observed near-ultraviolet spectral slopes of quasars in our 
sample lie in a quite wide range from —2.9 for the bluest one 
to for the reddest one, with a median value of —1.7. We plot 
FmirZ-Fnuv versas oi in Fig.[T] There is a fairly strong correlation 
between -Fmir /Fnuv and a with a Spearman ranking correlation 
coefficient Rs = 0.41, corresponding to a chance probability of 
Pr < 10~®. We fit this relation with a linear functior|j 

l0g(FMIR/FNUv) = fc ■ Q + 6, (2) 

where k and b are free parameters. For all 12483 quasars including 
in our primary sample, the best fit yields k — 0.24 ± 0.02 and 
6 = 0.45 ± 0.05. 

This correlation can be naturally inter preted as wave l ength- 
dependent dust extinction. For example, iRichards et al.l l l200lh 
showed that SDSS colors of quasars can be well represented by 
power laws of spectral indexes a = —1.5 ± 0.65 with a red tail 
possibly due to dust extinctions. Because extinction is larger at 
shorter wavelengths, it makes the observed optical-UV spectrum 
flatter (redder), and simultaneously raises observed Fmir/Fnuv. 
The large intrinsic scatters of a and Fmir /Fnuv dominate the 
variations at small a and small Fmir /Fnuv, thus make the cor- 
relation weaker for blue sources as seen in Figure [T] Given that 

^ see |http://www.sdss .org/dr7/products/spectra7spectrophotometry.html 
^ We use the IDL procedure robustjinefit.pro for ail the linear fitting in 
tliis work. 
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spectra slope a 

Figure 1. Correlation between Fmir/Fnuv and spectra slope a. The 
black dots represent the 12483 objects in our primary sample and red dots 
represent the 3792 objects with only upper limits of W3 or/and W4 magni- 
tudes. The green line displays the tendency predicted by SMC-like extinc- 
tion, only for those a > —1.7. The green dash line denotes the mean value 
of Fmir /Fnuv for quasars with a < —1.7. 



most blue quasars are likely not dust reddened, we check how the 
parameter k changes with different a ranges. For the 7700 sources 
with a > -1.7, we obtain ki.y = 0.29 ± 0.02, and for the 4505 
sources with a > —1.5, the result is fci.s = 0.33 ± 0.02. 

By comparing the behaviours of radio-loud and radio-quiet 
quasars in our sample, we find that this Fmir/Fnuv-Q relation is 
independent of radio activities. Also, for the 3792 sources that are 
removed from our primary sample, we calculate their Fmir using 
the same approach described in 12. 21 with the upper limits of W3 
& W4 magnitudes provided in the WISE catalog, and find that they 
share the same behaviour as for the primary sample (the red dots in 
Fig. [Hi. Therefore, we will not distinguish these different categories 
of quasars when analyzing extinction effect below. 

Previous studies suggested that dust extinction to quasars 
can be characterized by an SMC-like extinction curve (e.g. 
[Richards et aLll2003h iHopkins et aklliooi : lOaskell et al.ll2004h . To 
illustrate the extinction effect on Fmi r / Fnuv- q relation quantita- 
tively, we apply SMC-like extinction l lPei 19921) of different values 
of E{B — V) to an intrinsic blue QSO continuum, which is de- 
scribed by an intrinsic spectral slope /3 (Fa oc A'') in the NUV 
range and a certain fraction of reprocessed infrared emission, to 
get a set of faked quasar continuua. We measure the NUV spectral 
slope (a) and Fnuv for each faked continuum in exact the same 
way as that described above for real QSOs. It turns out that the re- 
lationship between a and E{B — V) can be well described by a 
linear function, 

E{B-V)^cia-P), (3) 

for E{B - V) < 0.3. Particularly, for P = -2.0, we have 
c = 0. 127. For such a.E{B — V) range, the dust extinction has little 
effect on the Fmir, which can thus be considered as a constant. The 
relation between log(FMiR/FNuv) and a for the mock quasars 
can be well described by a linear function with a slope of 0.322, in- 
dependent of the value of /3 for a wide range of —2.2 < /3 < —1.8. 
This slope is also very close to the observed one fci.y = 0.29 or 
fci,5 = 0.33. Therefore, our analysis supports that the observed 
Fmir/Fnuv-« relation is statistically attributed to the dust extinc- 
tion. And the scattering of this conelation is dominated by intrinsic 
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dispersions of the continuum slope and -Fmir/^nuv (green lines 
in Fig.[T](- 



3.2 Extinction Correction for -Fnuv 

In this subsection, we will correct Fnuv for dust extinction in a 
statistical way using the relation derived in the last section. First, 
we estimate the intrinsic distribution of log(FMiR/-FNUv) using 
blue quasars with a < —1.7. As mentioned above, in this regime, 
extinction is small, and variations of log(FMiR/-FNUv) and a are 
dominated by intrinsic scatters, rather than reddening effect, albeit 
the latter must be present. The distribution of log(-FMiR/JNUv) 
for blue quasars can be well described by a Gaussian function with 
a mean value of 0.010 (-Fmir/^nuv = 1-02), and a dispersion of 
a = 0.144. 

As we have shown in the last subsection, the dependence of 
Fmir I fj-iuv on observed spectral slope a can be explained by dust 
extinction. Therefore, it is possible to correct Fnuv in a statistical 
way for the extinction using the derived -Fmir/Fnuv-Q relation, 
with fitted spectral slope a as an extinction indicator. Assuming the 
mean Fmir/-Fnuv at different a is the same (1.02) as that for blue 
quasars, and following the FmirZ-Fnuv-Q relation derived in the 
last subsection, we calculate the correction factor for Fnuv: 



C 



3.5 X 10" 
1 



if Of > 
if a < 



-1.7 
-1.7. 



(4) 



In the following analysis, we still use the notation Fnuv to denote 
the NUV flux, while it has already been corrected for extinction 
effect. It should be pointed out that this correction is only in sta- 
tistical sense for the entire sample rather than for individual object 
since the intrinsic UV continuum slope /3 and Fmir/Fnuv have 
fairly broad distributions. Nevertheless, by applying different cut 
on a as the standard for blue quasars, we estimate the uncertainty 
of this correction for individual object is less than 10%. 



4 WARM DUST COVERING FACTOR AND ITS 
RELATION WITH AGN PROPERTIES 

4.1 Covering Factor of Warm Dust 

The observed MIR flux Fmir comes from three components: emis- 
sion of the warm and hot dust heated by AGN, young-star emis- 
sion at MIR band, and the emission of the accretion disc extending 
to MIR bands. Star-formation galaxies are characterized by strong 
PAH emission in MIR. Previous infrared spectroscopic studies sug- 
gest that PAH features are very weak in quasar spectra. Therefore, 
this component is not statistically important in our MIR range of 
rest-frame 3-10 /xm. We will neglect it in the analysis hereafter. 
We estimate the accretion disc component by extrapolation of the 
optical/UV power-law corrected for reddening effect. This could be 
considered as an upper limit considering spectral flattening in the 
near-infrared due to starlight contribution in the SDSS spectrunfl 
Numerically, the disc contribution to Fmir is 5.2% of Fnuv for 
/3 = -2.2, 9.3% of Fnuv for /? = -2.0 and 16.6% of Fnuv for 
/3 = —1.8. Since we focus on the MIR emission of the dust tori, we 
will subtract 0.093Fnuv from Fmir to remove the contribution of 
the accretion disc and the slight starlight. The remained value, also 



* Stai'light contribution in our co ntinuum windows is general small for lu- 
minous quasars in this paper (e.g. lStem & Laoil2012b . 



denoting by Fmir, should be considered as the MIR emission of 
the warm dust. 

Assuming that the NUV and MIR are isotropic, we can cal- 
culate the NUV luminosity Fnuv and MIR luminosity I/mir from 
Fnuv and Fmir, respectively. Next, we consider the bolometric 
coiTection (BC) for the total luminosity from optical to X-ray. 
Since we have already corrected for extinct ion effect, we will use 
the mean SED of optically blue quasars in [Richards et al. I ( l2006l) 
to calculate this correction. We calculate the Fnuv using the same 
method described in ij 12.21 and integrate log(;^/^) over log v from 
1 pxcL to 10 keV to obtain the total luminosity in this rang^. This 
yields BC — 4.34. Thus bolometric luminosities can be calculated 
by iboi ~ BC X Fnuv- Finally, we define the covering factor 
of warm dust as CFwd = L MiR/Lboi. Our de finition of covering 
factor is similar to that of Calderone et al. I (l2012l) . but different from 
theirs in the sense that we consider the IR flux in a fixed MIR range 
dominated by warm dust rather than the whole IR bands covered by 
WISE. Besides, we do not take into account anisotropic emission 
intentionally by this stage, while we will examine orientation effect 
in the last section. 



4.2 Dependence of CFwd on AGN Properties 

In this subsection, we consider the relations between CFwd and 
AGN properties: bolometric luminosity Lboi, black hole mass 
A/bh, and Eddington ratio L/ L-Edd- We present these relations for 
our sampl e in Fig. [2] We ad opt the virial black hole masses in the 
catalog of IShen et al.l ( 1201 ih with an additional extinction correc- 
tion. These masses are estimated using the width of Mgll A2798 
and U V continuum luminosity at 2800 A (Eqn. (2) in lshen et al.l 
l20Tlh . 



log( 



AfBH.vir 



a + b log ( 
+2 log( 



AF> 



1044 erg s 
FWHM, 



(5) 



kms~i 

with a = 0.740 and h = 0.62. IShen eTal] ( 1201 ih did not consider 
intrinsic extinction on the UV luminosity. In the following analysis, 
we will apply a correction by adding 

A = 1.588 * E{B - V) (6) 

to the value provided bv lShen et al ] (l20r l|). This correction is calcu- 
lated according to SMC reddening curve at 2800 A. And E{B — V) 
can be calculated from spectra slope (a) using Eqn. ([3) in § 13.11 
Bolometric luminosities and black hole masses for the entire sam- 
ple are also provided in the online table. 

First, CFwD is anti-correlated with Lboi. The Spearman rank 
correlation coefficient (Fs) of this correlation is —0.36, indicating 
the anti-correlation is very significant {Pr < 10^''). A linear fit 
yields 



log CFwD = (-0.23 ± 0.02) log Lboi/(erg ■ s" 
+ (9.79 ±0.04). 



(7) 



There is a significant anti-correlation between CFwd and 
black hole mass A/bh with a Spearman rank correlation coefficient 



^ We discount the IR luminosity here since most mi d- and far-infrare d 
emission is the re-radiation of the absorbed continuum dShang et al.ll201lh . 
But this still includes a small contribution from emission lines, which is 
doubly counted. 
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(Rs = -0.27, Pr < 10~^). A linear fit yields 

log CFwD = (-0.13 ± 0.02) log(A/BH/Mo) 

+(0.50 ±0.02). (8) 

Finally, no apparent correlation is found between CFwd and Ed- 
dington ratio, 

-±- = ^ (9) 

LEdd 1.5 X 1038 Mbh ^ ' 

{p = -0.008 and Pr = 0.34). 

We compare radio-quiet and radio-loud quasars on these three 
diagrams, and find that they are indistinguishable. Particularly, we 
have fitted related parametres separately for the radio-loud sub- 
sample. First, the median value of CFwd is consistent with that 
of the radio-quiet subsample within 10%. Second, the slopes of 
CFwD-iboi and CFwd-A/bh relations are —0.21 ± 0.02 and 
—0.15 ± 0.02, respectively. The distribution of CFwd, and these 
fitting slopes are consistent with their I a uncertainties. Therefore, 
we conclude that none of these three relations are affected by radio 
activities. 

Next, we consider the 3792 sources that only have upper lim- 
its of W3 and/or W4 magnitudes. From the left panel of Fig|2] it 
is obvious that these objects have systematically lower Lboi, while 
at the same time their CFwd (upper limits) are similar to those 
of quasars in our primary sample, so they shift to the left on the 
CFwD-iboi diagram. In other words, they show systematically 
lower CFwd, comparing to the quasars of similar luminosities in 
the primary sample. Because the fraction of these sources is small, 
we get almost an identical slope for CFwD-iboi relation for the 
entire sample. From the middle and right panels of Fig. |2] these 
quasars have similar distribution of black hole masses, but system- 
atically lower L/Lfidd than the primary sample. However, they fol- 
low the same CFwd-Mbh and CFwD-i/^Edd relations as the 
primary sample. 



5 DISCUSSION AND CONCLUSION 

5.1 On Radio-loud Quasars 

Approximately 10% objects in our sample are radio-loud quasars. 
As we have discussed earlier, they exhibit the same correlations 
as radio-quiet quasars. Therefore, the majority of the radio-loud 
quasars have the very simil ar structure of accretio n di sc and dusty 
torus. This is consistent with lRichards et aLl ( l2006h and lShang et al.l 
( |20T l|) that mean SEDs of radio-quiet and radio-loud quasars in 
NUV and MIR are quite similar. However, it is of necessity to point 
out that a small fraction of the radio-loud quasars within our sam- 
ple are blazars, in which the relativistic jets are beamed toward 
us. In this case, the relativistic jets would contribute significantly 
or even dominate the observed flux at near-ultraviolet and mid- 
infrared bands. Therefore, NUV and MIR emission in these ob- 
jects will no longer represent the primary emission of the accretion 
discs and re-emission of the dusty tori. Thus our definition of warm 
dust covering factor is no long appropriate for these blazars. Nev- 
ertheless, the fraction of blazars in our whole sample is q uite small 
(~ 1 % of type I quasars according to the argument of |Padovar3 
l2007h . they will definitely not affect our conclusion. 

5.2 Comparison with Previous Studies 

In this work, we use the mid-infrared to near-UV luminosity ratio 
to estimate the covering factor of warm dust, and find a decrease of 



CFwd with increasing Lboi. Our results are in qualitative consis- 
tenc y with previous works. 

iMaiolino et al. I l l2007l) estimated the dust covering factor for 
a few tens of QSOs using monochromatic flux ratio of 6.7/im to 
5100A, which is similar to our method in the spirit. However, their 
sample was constructed by merging two groups of QSOs from to- 
tally separate redshift ranges. Moreover, their analogous CF-Lboi 
correlation is far less significant when only considering either one 
group of objects with close redshifts. Therefore, their result might 
be hampe red by possible evolu tion effect. A similar trend was ob- 
tained bv lTreister et al.l l l2008h . who made use of nearly one hun- 
dred objects constructed from different surveys. In comparison with 
these works, our sample is much larger, and is uniformly selected 

in a relatively small redshift range^ 

More recently, [Mor & Tra khtenbroj (|20T ij) used the prelim- 
inary data release of WISE to study the hot dust component for 
a large and uniformly selected sample. Their definition of hot 
dust covering factor is quite similar to ours, and their results also 
revealed a dependence of hot dust CF on bolometric luminos- 
ity. Their sample covers a very wide redshift range and their hot 
dust luminosities were estimated using a certain dust model. In- 
stead , we take a model-inde pendent approach in our work. Simi- 
larly, |Calderone|et^](|20l2) studied the dust covering factor for all 
radio-quiet AGNs constructed from WISE and SDSS within a small 
redshift range. These authors reported a similar trend of CF-Lboi 
relation by dividing the whole sample into three bins. Comparing 
with these studies, our work has taken the extinction effect into 
account. Also, we used the extinction-corrected broad band NUV 
luminosity, which is a better estimation of bolometric luminosity. 
Furthermore, we find that CF is significantly correlated with black 
hole mass, but is independent of Ed dington ratio . The latter results 
are consistent with those reached bv lCaol l l2005h for PG quasars. 

In the literature, many authors used the fraction of type 2 ob- 
jects as an indicator of total d ust covering facto r and reached a sim- 
ilar co nclusion. For instance, ISimpsonl fcOOSh and iHao & Strauss! 
( I2OO4I) found a similar dependence of type 2 fraction on the 
[O III] luminosity using a magnitude-limited sample from the 
SDSS spectroscopic survey. Similar results were reported by other 
authors using sample co nstructed from flux-limited X-ray sur- 
vey (e.g. Hasinge3l2008l) or a combination of multi-band survey 



jHatziminaoglou et alj|2009l) . However, some o t her au thors have 
reported the opposite results. [Lawrence & Elvis] (I2OIOI) suggested 
that there is no correlation between type 2 fraction and [O III] lu- 
minosity after re moving LINER s from the definition of type 2 ob- 
jects, and so did lLu et"^ 1 I2OIOI) by taking into account of various 
selection effects. A caution should be taken when comparing our 
results with those derived from the fraction of type 2 AGNs. We 
only consider the covering factor of warm dust, which is merely 
a portion of the obscuring material. Furthermore, a fraction of the 
type 2 objects are ob scured by dust in their host gal axies rather than 
dust tori themselves jGoulding & Alexandeill2009l) . 

5.3 Effects of Anisotropic Continuum 

It should be emphasized that up to now, our whole analysis is based 
on the assumption that the NUV and MIR emission of quasars are 
isotropic. However, the NUV emission from a geometrically thin 



lar directions than at equatorial directions (e.g. 


Shakura & SunvaevI 


I1973I; iLaor & Netzeilll989l iFukue & Akizukil 


20061). If this is the 



case, our estimation of bolometric l uminosity would be biased de- 
pending on the viewing angle (cf . iNemmen & BrothertonI l2010t 
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Figure 2. Relations between warni dust covering factor CF^D and AGN properties: bolometric luminosity Lbol (l^ft panel), black hole mass A/bh (middle 
panel), and accretion rate ibol/^Edd (right panel). The meanings of black and red dots are exactly the same as in figure 1. The blue dots represent radio-loud 
quasars in the primary sample. The green dash lines represent the hnear fit for these relations. 



iRunnoe et al and our definition of covering factor would 

certainly be inaccurate. For example, if an object is viewing rel- 
atively edge-on, we will underestimate its bolometric luminosity 
Lboi while overestimate its covering factor CFwd. In a statisti- 
cal view, this can result in an anti-coiTelation between CFwd and 
Lboi, the same tendency as we have discovered from our sample. 

To illustrate this more clearly, we construct a simple model 
to simulate the effect of different viewing angles. The simulation 
involves the same number of faked quasars as for our sample. 
We adopt three assumptions: (1) The intrinsic CFwd and Lboi 
are independent, and obey log-normal distributions respectively. 
(2) The inclinations are random in a range of 1/2 < fi < 1, 
where n = cos 9 (the viewing angle). A lower limit on fi is set, 
since we require the line-of-sight must lie in the open angle of 
the torus, otherwise the NUV emission would be totally absorbed. 
The value 1/2 is corresponding to a type 2 to type 1 ratio of 1 : 10 
jLawrence & Elvisl2010l : lReves et al.l2008h . (3) The observed iboi 
and CFwd for each object are calculated from observed NUV flux 
and NIR/NUV flux ratio viewing from a certain direction. The de- 
pendence of NUV flux on viewing angle used in our model is 



/(m) = /o X -), 



(10) 



following the results of iFukue & Akizukil ( |2006|) who have taken 
into account limb darkening effect. Moreover, we require that the 
observed Lboi and CFwd have exactly the same distribution as for 
our sample. 

We plot our simulated CFwD-Lboi relation in Fig. [5] in com- 
parison with the same relation of the sample. It turns out that 
our simulation obtains a similar but more compact relation. This 
seems to indicate that the observed CFwD-Lboi may result from 
anisotropic NUV continuum. 

However, we stress that our simulation is initial and needs to 
be refined in several aspects. First, the assumption that the viewing 
angles are uniformly distributed is not satisfied by a flux-limited 
sample, such as the SDSS quasar sample used here. For example, a 
quasar close to the flux limit may move out of this limit if it is ob- 
served at a higher inclination. Thus it will be missed in the sample. 

^ We set the lowest ^ to the value for the average torus oppening angle, 
rather than that giving by the warm dust covering factor of each object be- 
cause the coveting factor of thick cold and hot dust is not known. 
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Figure 3. Above: the CFwD-Lboi relation for the piimary sample. Bot- 
tom: the same relation for the simulated sample. 



In other words, the sample is biased against least luminous quasars 
at large inclinations. We note that the continuum flux decreases by 
a factor of 2.9 from fi — 1 to /i — 1/2, so this selection effect 
will significantly affect the inclination distribution of quasars with 
luminosities of a factor of 2.9 above the minimum value. The net ef- 
fect will weaken the CFwD-Lboi correlation. In principal, we can 
incorporate this selection effect into the simulation by using the 
quasar luminosity function. However, we should understand that 
dust extinction is also coupled with the selection effect in the real 
case. Furthermore, dust extinction may or may not correlate with 
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the inclination. If the absorber is the continuous extension of the 
dusty torus, there will certainly be a correlation; while on the other 
hand, if the absorption is due to the interstellar medium of the host 
galaxy, there may be no correlation. Because of lack of such de- 
tailed knowledge, we did not carry more complicated simulations. 

Second, IR emission is also likely to be anisotropic due 
to the radiation transf er effect. According to the calculations of 
iNenkova et al.l ( l2008al lbl). for a typical clumpy dusty torus model, 
the MIR flux is ~ 50% weaker viewing from a 45 degree inclina- 
tion than viewing from face-on direction. This factor is very simi- 
lar to the anisotropy of NUV flux. In that case, the CFwD-^boi 
relation will be much shallower due to the decrease of infrared 
flux viewing from higher inclinations. Hence the correlation will 
be weakened. 

Therefore, although our simulation suggests that the observed 
CFwD-iboi relation may be attributed to the inclination effect, 
there is fairly possible that the relation is intrinsic. A detailed study 
combining anisotropic radiation from both the discs and dust tori 
is very favourable in the future. Also, continuum variation on time 
scales of years will cause similar trend due to a delayed response 
of dust emission. The fact that CF is con'elated with black hole 
mass but not correlated with Eddington ratio further indicates that 
luminosity dependence of CF is not solely due to the anisotropic 
emission alone. 

To summarize, we use a complete sample in a narrow redshift 
range to study the warm dust emission of ~ 16, 000 quasars. We 
correct for the dust extinction effect in a statistical way for the first 
time in this kind of study. We find that CFwd is con'elated with 
Z/boi and AIbh, while not with Eddington ratio. Monte Carlo sim- 
ulation suggests that CFwD-iboi relation may be significantly af- 
fected by anisotropic disc emission. But it may not be the dominant 
factor for this relation. 
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APPENDIX 

We present the primary data of all 16275 sources in our sample 
with an ASCII table as the supplementary material of this paper. 
The entire table is available online. Here we offer a small portion 
of the table for a guidance (Table [T}. 
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Table 1. Primaiy data of our sample. 



Name 


Redshift 


Spectra slope 


NUV flux 


MIR flux 


CFwD 


Bolometric luminosity 


Black hole mass 


Flag'' 


FIRST 


(SDSS J) 


(z) 


(-a) 


(i^NUv") 


(-Fmir") 




log(Lboi/(ergs~^)) 


log(MBH/AfQ) 




flag'^ 


000013.14+141034.6 


0.958 


1.657 


2.773e-13 


2.617e-13 


0.217 


45.757 


8.800 





_j 


000024.02+152005.4 


0.989 


1.744 


2.820e-13 


2.044e-13 


0.167 


45.798 


9.181 


1 




000025.21+291609.2 


0.924 


1.718 


8.433e-14 


9.096e-14 


0.249 


45.201 


8.216 







000025.93+242417.5 


1.156 


1.296 


1.487e-13 


2.077e-13 


0.322 


45.688 


9.310 


1 




000026.29+134604.6 


0.768 


1.691 


3.206e-13 


1.291e-13 


0.093 


45.582 


8.670 


1 




000028.82-102755.7 


1.152 


1.277 


4.372e-13 


2.150e-13 


0.113 


46.152 


9.449 







000031.86+010305.2 


1.093 


1.086 


2.274e-13 


3.818e-13 


0.387 


45.812 


8.552 








000042.89+005539.5 


0.951 


1.805 


9.408e-13 


1.033e-12 


0.253 


46.279 


8.844 








000057.67-085617.0 


1.095 


1.460 


8.727e-13 


8.643e-13 


0.228 


46.397 


8.778 








000123.96-102458.1 


1.058 


1.616 


6.734e-13 


6.254e-13 


0.214 


46.248 


9.338 








000123.98+284249.4 


0.936 


0.269 


2.419e-13 


2.935e-13 


0.280 


45.672 


9.255 





-1 


000140.70+260425.5 


0.765 


1.385 


3.656e-13 


3.700e-13 


0.233 


45.636 


8.670 





-1 



The full table including all 16275 sources is available onHne as supplementary material. A portion is shown here for guidance of its content. 
" The flux is in unit of erg cm~^ A~^. 

Here -Fnuv refers to the value after extinction correction (see §|3) and Fmir refers to the value after discounting the disc contribution (see ^ 14. It . 
Flag: 0=12483 source in our primary sample; l=the rest 3792 sources which only have upper limits for W3 or W4 magnitudes. 

First flag: the same as column 15 in the catalog of^hene£al^^201^. -l=not in FIRST print; 0=FIRST undetected; l=core dominant; 2=lobe dominant. 
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